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We predict that the recently discovered quasi-one dimensional superconductors, A 2 CrsAs 3 (A=K,Rb), pos¬ 
sess strong frustrated magnetic fluctuations and are nearby a novel in-out co-planar magnetic ground state. The 
frustrated magnetism is very sensitive to c-axis lattice constant and can thus be suppressed by increasing pres¬ 
sure. Our results qualitatively explain strong non-Fermi liquid behaviors observed in the normal state of the 
superconductors as the intertwining between the magnetism and superconductivity can create a large quantum 
critical region in quasi-one dimensional systems and also suggest that the materials share similar phase diagrams 
and superconducting mechanism with other unconventional superconductors, such as cuprates and iron-based 
superconductors. 


PACS numbers: 74.70.-b, 74.25.Ha, 74.20.Pq, 74.20.Rp 

One of major challenges in condensed matter physics is to 
understand the role of electron-electron correlation in uncon¬ 
ventional superconductors. The effect of electron-electron in¬ 
teraction becomes more important as the dimension of a sys¬ 
tem is lowered. Indeed, many unconventional superconduc¬ 
tors discovered in the past are quasi-two dimensional(Q2D) 
electron systems. The superconductivity in these unconven¬ 
tional superconductors appears in a vicinity to a magneti¬ 
cally ordered state. Magnetic fluctuations which are caused 
by electron-electron interaction have been widely considered 
to be responsible for superconductivity and many non-Fermi 
liquid behaviors in normal states. 

While there are many representatives of Q2D unconven¬ 
tional superconductors, it has been difficult to find one 
in quasi-one dimensional(QlD) systems even if the effect 
of the electron-electron correlation is expected to be en¬ 
hanced further. The Q1D superconductors discovered pre¬ 
viously, including Bechgaard salts[l, 2] , Tl 2 MogSe6[3] and 
Lio.gMogOi7[4—7], are not attributed to 3d-orbital electrons 
which can exhibt strong electron-electron interaction. 

Very recently, two novel Q1D materials K 2 Cr:>As:j[8] and 
Rb 2 Cr 3 As3[9] have been synthesized and found to be super¬ 
conducting below the transition temperature 6.1 K and 4.8 K 
respectively. The structure of A 2 Cr 3 As3 (A=K,Rb) is charac¬ 
terized by one-dimensional (Cr 3 As 3 ) chains (Fig. 1(a)), which 
contain Cr fi distorted octahedral clusters. The alkali metal 
ions are intercalated between the (CT3AS3) chains. Both new 
materials show strong non-fermi liquid behaviors in normal 
states, as well as unconventional superconducting properties 
in superconducting (SC) states. Moreover, just like cuprates 
and iron-based superconductors, the electronic physics in 
these new materials are likely attributed to 3d-oribtals of Cr 
atoms. Therefore the material may exhibit strong magnetism 
and electron-electron interaction. 

In this paper, we show that the new materials can be a crit¬ 
ical representative of Q1D unconventional superconductors 
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where the superconductivity emerges in a vicinity to a novel 
magnetically ordered state. We predict that the materials are 
characterized by strong frustrated magnetic fluctuations and 
are nearby a novel in-out co-planar(IOP) magnetic ground 
state. The magnetism can be described by a minimum effec¬ 
tive model with three magnetic exchange parameters: the anti¬ 
ferromagnetic ./| and J\ between two nearest-neighbor (NN) 
Cr atoms and the ferromagnetic J 2 between two next NN Cr 
atoms along c-axis. The frustrated magnetism is very sensi¬ 
tive to c-axis lattice constant and can thus be suppressed by in¬ 
creasing pressure. The results suggest that the materials host a 
typical phase diagram similar to those of the Q2D unconven¬ 
tional superconductors, such as cuprates [10] and iron-based 
superconductors [11]. The new materials can be ideal systems 
to study the intimate relations between magnetism and super¬ 
conductivity since a Q1D model can be solved theoretically 
with high controllability. 

It is known that calculations based on the Density Func¬ 
tional Theory (DFT) can successfully identify the possible 
magnetic ground state in some correlated electron systems. 
Typical examples are the iron-based superconductors where 
theoretical calculations consistently agree well with experi¬ 
mental measurements[12-18] even if the magnetic origin is 
still debatable [19]. Here, we deploy similar DFT calculations. 
Our DFT calculations employ the projector augmented wave 
(PAW) method encoded in Vienna ab initio simulation pack- 
age(VASP) [20-22], and generalized-gradient approximation 
(GGA)[23] for the exchange correlation functional is used. 
Throughout this work, the cutoff energy of 450 eV is taken 
for expanding the wave functions into plane-wave basis. The 
number of these k points are 6 x 6 x 13 for K 2 Cr 3 As3 and 
Rb 2 Cr 3 As 3 [24]. The GGA plus on-site repulsion U method 
(GGA+U) in the formulation of Dudarev et a/. [25] is em¬ 
ployed to describe the electron correlation effect associated 
with the Cr 3 d states by an effective parameter t/ e //. The val¬ 
ues of (7=2.3 eV and ,7=0.96 eV on Cr are adopted for our 
GGA+U calculations[26]. We relax the lattice constants and 
internal atomic positions, where the plane wave cutoff energy 
is 600 eV. Forces are minimized to less than 0.01 eV/A in the 
structural relaxation. 
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FIG. 1: Schematic view of the structure of A 2 Cr 3 As 2 (a), (b) and 
(c) show the all-in and in-out noncollinear magnetic states. The ex¬ 
change coupling parameters Ji, J[ and J 2 are defined in (d). The 
orange and purple spheres represent Crl and Cr2. 


TABLE I: Experimental and optimized structural parameters of 
K 2 Q" 3 As 3 using GGA in the paramagnetic phase. Deviations be¬ 
tween the optimized and experimental values are given in parenthe¬ 
ses in % . 



GGA 

EXP 

a( A) 

10.113(+1.3) 

9.983 

c(A) 

4.147(-1.98) 

4.230 

Crl-Crl/Cr2-Cr2(A) 

2.498; 2.588 

2.615; 2.691 

Cr 1 -As l/Cr2-As2(A) 

2.513; 2.494 

2.51; 2.49 

Cr 1 - As2/Cr2-As 1 (A) 

2.522; 2.506 

2.516; 2.506 

ai /a 2 (°) 

59.6; 62.51 

62.8;65.4 

/W 2 O 

60.8; 60.87 

62.8;62.9 

7i/72 (°) 

110.6; 111.6 

114.4; 115.2 


Due to the asymmetric distribution of alkali metal ions and 
the absence of inversion symmetry, there are two kinds of 
Cr ions in K 2 Cr 3 As 3 (Rb 2 Cr 3 As 3 ): Crl and Cr2. The Crl 
ions are surrounded by six inplane A ions and the inplane 
bondlength is 2.615(2. 57)A. The Cr2 ions are surrounded by 
three inplane alkali ions and the inplane bondlength is a bit 
longer, 2.69(3.16)A. From the number of A ions surrounding 
Cr ions, it is expected that the CrlAs layer will obtain more 
electrons than the Cr2As layer. 

The optimized and experimental structural parameters of 
K 2 Cr 3 As 3 are summarized in Table I. We find that the lattice 
constants are comparable with experimental values but the Cr- 
Cr bond lengths are underestimated by 0.1 A. Thus, the Crl¬ 
As 1 -Crl(ai)/Cr2-As2-Cr2(a 2 ), Crl - As 1 -Cr2(/3-| )/Cr2-As2- 
Crl(/3 2 ) and Crl-As2-Crl(7!)/Cr2-Asl-Cr2(7 2 ) angles are 
also underestimated. Similar cases have been also noted in 


the studies of iron based superconductors [17, 27, 28] and it 
may be related to the strong spin fluctuation, which is beyond 
DFT calculation. We adopt the experimental parameters in the 
following calculation unless otherwise specified. 

The band structure, the density of states(DOS) and Fermi 
surfaces for K 2 Cr 3 As 3 have been calculated in Ref.29. Our 
results as shown in Appendix A are similar to their calcula¬ 
tions. In general, the 3d states of Cr are located from -2.5 eV 
to 2.0 eV and the states near the Fermi level are mainly at¬ 
tributed to Cr d z 2 , d xy and d x 2_ y 2 orbitals. The As 4 p states 
mainly lie 1.0 eV below the Fermi level and hybridize strongly 
with the Cr 3 d states. Along k z direction, the bands are very 
dispersive due to the strong coupling between d oribials along 
z direction. The bands attributed to Cr2 d xy and d x 2_ y 2 or¬ 
bitals are very close to half-filling. As the separation between 
chains is large, the inplane band dispersion is very small but 
not negligible. However, we want to point out an important 
feature that is ignored in Ref.29. Crl ions have more d elec¬ 
trons than Cr2 ions, which is consistent with analysis from the 
crystal structure. We will show later that this difference also 
results in different magnetic moments at Crl and Cr2 sites in 
magnetic states. 

The band structure of Rb 2 Cr 3 As 3 , shown also in Appendix 
A, shares many common features with that of K 2 Cr 3 As 3 , but 
they also differ in some details. As the radius of Rb is big¬ 
ger than K, CrlAs layer can get more electrons and d xy and 
d x 2_ y 2 bands of Crl are fully occupied. Meanwhile, the d xy 
and d x 2_ y 2 bands of Cr2 are much less occupied. It leads to 
great difference between the 3D Fermi surfaces of the two 
materials. Another difference is that there is an additional 
electron Fermi surface near A(0, 0,7r) point in Rb 2 Cr 3 As 3 , 
which is attributed to d xy and d x 2_ y 2 orbitals of Cr2. Finally 
the d xz and d yz orbitals are much close to the Fermi level in 
Rb 2 Cr 3 As 3 . The calculated N (Ep) in K 2 Cr 3 As 3 (Rb 2 Cr 3 As 3 ) 
is 8.76(9.13) eV -1 /f.u.. The calculated Pauli susceptibility 
and specific heat coefficient are \o = 2.83(2.95) x 10 -4 
emu/mol and 7 = 20.7(21.5) mJ/(A' 2 mol). The calculated 7 
is about only one third of the experimental value in K 2 Cr 3 As 3 , 
suggesting strong correlation in these systems. 

Now we focus on the magnetic properties. We consider 
four possible collinear magnetic states, the paramagnetic state 
(PM), the ferromagnetic (FM) state, the interlayer antiferro¬ 
magnetic (AFM) state and the up-up-down-down (tt4A) mag¬ 
netic state. Due to strong magnetic frustration, we also con¬ 
sider two additional co-planar antiferromagnetic states: all-in 
magnetic state (Fig. 1(b)) and in-out magnetic state (Fig. 1(c)). 
We perform calculations with spin orbital coupling (SOC) and 
the calculated magnetic moments and the total relative ener¬ 
gies of above magnetic states are summarized in Table.II. The 
IOP state is the ground state in K 2 Cr 2 As 2 (Rb 2 Cr 2 As 2 ), with 
a large energy gain of 39(762) meV/cell relative to the PM 
state. The magnetic moments are 0.90(1.75) and 0.94(2.34) 
ji [> on Crl and Cr2 sites, respectively. The initial FM state 
converges to a PM or AFM state, indicating that the inter¬ 
layer magnetic coupling J[ as indicated in Fig. 1(d) is antifer¬ 
romagnetic and relatively strong. In the AFM state, the energy 
gain is 25(486) meV/cell and the three magnetic moments at 
Crl or Cr2 site are different. It manifests that the intralayer 










3 


magnetic couplings among three Cr atoms J-\ is AFM. Thus, 
strong magnetic frustration exists in A 2 Q 3 AS 3 . Our results 
are different from those in Ref. [29] which missed to identify 
the IOP as the true magnetic ground state in their calculations 
To confirm our results, we have performed all-electron calcu¬ 
lations and find that relative energies and magnetic moments 
for collinear states are rather close to those in Table II. 

The IOP state gains energy rapidly and becomes very robust 
if we include additional U, the onsite electron-electon corre¬ 
lation. The results from GGA+U calculations are given in 
Table. III. The IOP state in K 2 Cr 2 As 2 (Rb 2 Cr 2 As 2 ) has a large 
energy gain of 392(466) meV/cell relative to the AFM state. 
The magnetic moments are greatly enhanced to 2.49(2.76) and 
2.56(2.92) /is on Crl and Cr2 sites as well. It is also impor¬ 
tant to note that the magnetism is always found to be stronger 
in Rb 2 Cr 3 As 3 than in K 2 Cr 3 As 3 . In the presence of U, as 
the calculation can be converged for different magnetic states, 
we can extract the magnetic exchange parameters from their 
energy differences [30]. In the GGA+U calculations, the local 
magnetic moments on Cr are very close to each other. Thus, 
we can can extract the magnetic exchange parameters from 
their energy differences within the Heisenberg Model. The en¬ 
ergies contributed by magnetic interactions in these five mag¬ 
netic states per cell are written as. 


Efm/S 2 
Eafm/S 2 

Enu/s 2 

Eall—in / S 2 

Ein—out/S 


6.7, + 12.7] + 6J 2; 

6 J 1 - 12 + 6 J 2 , 

6 J 1 - 6 J 2 , 

—3«7i + 6 J'i + 6J 2) 
—3Ji — 6J[ + 6 J 2 . 


(1) 


Using the relative energies and the above equations, we can 
obtain four equations. But there are only three exchange pa¬ 
rameters. For this overdetermined system of equations, we 
can obtain these parameters by using a least-squares tech¬ 
nique. The estimated exchange couplings are Ji = 0.090 
eV/S 2 , J[ = 0.060 eV/S 2 , J 2 = -0.010 eV/S 2 forK 2 Cr 3 As 3 
and Ji = 0.075 eVAS 2 , J[ = 0.045 eV/S 2 , J 2 = -0.020 
eV/.S ' 2 for Rb 2 Cr 3 As 3 in the GGA+U calculations. The NN 
exchange couplings Ji and ,7| are strongly AFM. However, 
the next NN exchange coupling J 2 is ferromagnetic. The IOP 
state saves energy from all these effective exchange couplings. 

Fig.2(a) and (b) show the calculated band structure and 
DOS in the IOP magnetic state of K 2 Cr 3 As 3 . The material 
remains metallic in the presence of the static magnetic order. 
The bands near the Fermi level split due to SOC and magnetic 
order but the orbital characters have little changes compared 
with those in PM state. This splitting depends on orbital char¬ 
acters of the bands: it is larger for d xy ,d x 2_ y 2 and d xz , d yz 
bands than for d z 2 bands. 

As the size of Rb atoms is much larger than the size 
of K atoms, the stronger magnetism in Rb 2 Cr 3 As 3 than in 
K 2 Cr 3 As 3 suggests that applying pressure can suppress mag¬ 
netic fluctuations. To simulate the effect of pressure, we in¬ 
vestigate the magnetism as a function of the lattice constant 
c. Fig. 3(a) shows the calculated magnetic moments and rel¬ 
ative energies of the IOP magnetic state(relative to PM state) 
as a function of c. As the decrease of c, which is equivalent 



FIG. 2: Band structure and DOS for K 2 CT 3 AS 3 in in-out co-planar 
(IOP) magnetic state. 


TABLE II: The total energies for different magnetic states of 
A 2 CT 3 AS 3 . They are given relative to the total energy of the PM 
state and the unit is meV/cell. The magnetic moments are given in 

fJ-B- 


A 2 & 3 AS 3 


relative energy (meV/cell) 

Mcr 1 

Mcr2 

A=K 

PM 

FM 

AFM 

tm 

all-in 

in-out 

0 

to PM 

-25 
to PM 
-3(to in-out) 

-39 

0 

0.43 

0.06 

0.90 

0 

0.53 

0.35 

0.94 


PM 

0 

0 

0 


FM 

-589(to AFM) 

0.40 

2.29 

A=Rb 

AFM 

-486 

1.26 

2.35 


tm 

-525 

0.01 

1.83 


all-in 

-594(to in-out) 

0.34 

2.15 


in-out 

-762 

1.75 

2.34 


to the increase of pressure, the magnetic moments decrease, 
as well as the energy gains. When c = 4.06 A, the energy 
gain reaches zero, indicating the vanish of IOP magnetic or¬ 
ders. We also directly calculate the pressure effect using the 
GGA+U calculation. Consistent results are obtained as shown 
in Fig.3(b). These results allow us to conjecture a phase dia¬ 
gram of these new superconductors similar to the one shared 
by many Q2D unconventional superconductors as sketched in 
Fig. 4(a). The non-Fermi liquid behaviors observed in experi¬ 
ments can be naturally explained as a result of strong magnetic 
fluctuations in the critical region. 

We can construct a minimum effective Hamiltonian to de¬ 
scribe the magnetism for a Q1D chain in A 2 C 13 AS 3 with the 
above mentioned magnetic exchange couplings. The Hamil¬ 
tonian is 

Em — J\ ^ ^ T <7 3 ^ ^ * S ^/3 

*.«/9 (ij),ci0 

+ J2 E Szet ' Sj/ 3 1 

{(ij)),ctP 


( 2 ) 
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TABLE III: The total energies for different magnetic states of 
AoCr 3 As 3 with GGA+U calculations. They are given relative to the 
total energy of the AFM state and the unit is meV/cell. The magnetic 
moments are given in ys- 


A 2 Cr 3 As3 


relative energy (meV/cell) 

Mcr 1 

Mcr2 


FM 

+ 1495 

2.47 

2.58 


AFM 

0 

2.48 

2.41 

A=K 

tm 

+854 

1.50 

2.44 


all-in 

+311 

2.35 

2.43 


in-out 

-392 

2.49 

2.56 


FM 

+ 1045 

2.09 

2.84 


AFM 

0 

2.70 

2.97 

A=Rb 

tm 

+792 

2.72 

2.90 


all-in 

+ 184 

2.90 

2.96 


in-out 

-466 

2.76 

2.92 




FIG. 3: The relative energies and magnetic moments as function of 
lattice constant c (a) and pressure (b). In the case of pressure, we 
performed structure relaxation and GGA+U calculations. 


where a,/? label sublattices and (ij) and ((ij)) denote the 
NN and next NN pairs. With Ji,,/[ being AFM and J 2 be¬ 
ing FM, Hm has the IOP magnetic ground state as a classical 
spin model or in the large S-limit. The magnetic excitations 
can be calculated by employing the Holstein-Primakoff trans¬ 
formation within the linear spin-wave approximation. As the 
ID chain described by Hm has a C' 3 ,, symmetry with six Cr 
atoms in the unit cell, there must be three acoustic and three 
optical spin excitations. Each of these three modes must be 
composed of one A\ and two degenerated E modes. The de¬ 
tailed analytic results of the spin waves are given in Appendix 
B. A typical spin wave spectrum is plotted in Fig. 4(b). The 
acoustic spin wave dispersion shows a linear dependence with 
small k z , which indicates the AFM nature of the magnetic ex¬ 
change couplings. 

Although there is no direct experimental measurement on 
magnetic properties in these materials, the magnetic exchange 
coupling parameters obtained from our calculation are con¬ 
sistent with recent experimental measurements in CrAs[31], 
a MnP-type orthorhombic crystal structure and also a super¬ 
conductor under pressure[32]. The magnetic order in CrAs is 
a double helimagnetic structure. Along the one-dimensional 
helimagnetic structure, the measurements have shown that it 



FIG. 4: (a) The conjectured phase diagram with pressure of 

A 2 Cr 3 As 3 . (b) The calculated spin-wave dispersion of A 2 Cr 3 As 3 
along k z direction, with Jj=0.7Ji and J2=-0.15Ji. 


is AFM between two NN Cr atoms that has a Cr-As-Cr angle 
~ 72° and FM between two next NN Cr atoms that has a Cr- 
As-Cr angle ~ 123°. These two magnetic exchange couplings 
in CrAs resemble Ji(J[) and J 2 . It is interesting to note that 
the exchange couplings between two Cr atoms are against the 
simple Goodenough-Kanamori-Anderson (GKA) rules which 
would suggest opposite signs for Ji(J[) and J 2 . Such a vi¬ 
olation suggests that multi-magnetic mechanisms may work 
together in CrAs based structures and J 2 may include signifi¬ 
cant contribution from the double-exchange mechanism since 
the average valence of Cr atoms in A 2 Cr 3 As 3 is 2.3. 

We have ignored the Dzyaloshinskii-Moriya(DM) term that 
can be induced by spin orbital couplings in the absence inver¬ 
sion center between two Cr atoms. In CrAs, the DM term 
is important in determining the incommensurate wavevector 
in the double helimagnetic structure. However, in A 2 CP 3 AS 3 , 
the effect of the DM is greatly reduced as the Cr atoms form a 
tight octahedral cluster. This is consistent with our calculation 
results that magnetic moments tend to be in the xy plane. 

Our results suggest that similar to cuprates and iron-based 
superconductors, the magnetic fluctuations are responsible for 
the superconductivity in A 2 Ci' 3 As 3 . If this is the case, it is 
interesting to ask what type of pairing is favored from the 
specific magnetic fluctuations predicted here. We argue that 
the favorite pairing is likely to be spin-triplet pairing from the 
magnetic fluctuations. As shown in Appendix A, the Crl and 
Cr2 atoms make major contributions to different bands near 
the Fermi level. Thus, as the pairing is expected to be dom¬ 
inated by the intra-band pairing, in the real space picture the 
majority pairing takes place within each sublattice that only 
contains one type of Cr atoms. The effective magnetic fluctu¬ 
ations within each sublattice along the the chain direction con¬ 
nected by J 2 are clearly FM from the above results. Therefore, 
if we ignore the parity breaking caused by alkaline atoms, the 
pairing from the magnetic fluctuations is very likely to be a 
p -wave spin triplet as the correlation effect forbids any onsite 
pairing. If we consider that the pairing is determined by the 
local FM exchange coupling[33], the induced gap function , 
A oc sink z , which should be characterized by line nodes on 
the Fermi surface in the k z = 0 plane in reciprocal space. 

In summary, we predict the Q1D superconductors, 
A 2 Cr 3 As 3 , are close to a novel in-out co-planar magnetic or¬ 
dered state and share a typical phase diagram similar to those 
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of Q2D high temperature superconductors, cuprates and iron- 
based superconductors. The prediction qualitatively explains 
the non Fermi-liquid behaviors observed in these systems and 
suggests that the superconductivity in these systems are driven 
by electron-electron correlation effects. We also predict that 
T c can be maximized in these systems by applying certain ex¬ 
ternal pressure. The new materials can be an ideal Q1D sys¬ 
tem to understand the intimate relation between magnetism 
and superconductivity. 
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Appendix A: Electronic structures of A 2 Cr 3 As 3 

The band structure and density of states(DOS) for 
K 2 Cr 3 As 3 with experimental structural parameters are shown 
in Fig. 5 (a) and (b),(c), respectively. From the band struc¬ 
ture and DOS, we find that Crl ions have more d electrons 
than Cr2 ions, which is consistent with analysis from the crys¬ 
tal structure. It will lead to different magnetic moments at 
Crl and Cr2 sites in magnetic states. The Fermi surfaces are 
shown in Fig.7(a), which are similar to those in Ref.29. The 
band structure and DOS for Rb 2 Cr 3 As 3 , shown in Fig. 6 , are 
similar to those of K 2 Cr 3 As 3 . The main difference is that the 
3D Fermi surface is quite different from that of K 2 Cr 3 As 3 and 
there is an additional electron Fermi surface near A( 0 , 0 , 7 t) 
point, as shown in Fig.7(b). 



FIG. 7: Fermi surfaces of K 2 G 3 AS 3 (a) and Rb 2 Cr 3 As 3 (b) with 
experimental parameters in the paramagnetic state. 


Appendix B: Spin excitations 

The model of Cr octahedral clusters is shown in Fig. 8 . The 
magnetic excitations in the effective magnetic Hamiltonian, 
Eq.2 , for A 2 CT 3 AS 3 can be calculated by using the Holstein- 



Primakoff transformation. 


7 IG. 5: 
>aramet 


Band structure and DOS of K 2 Cr 3 As 3 with experimental 
rs in the paramagnetic state. 



^ai — y 2jS 0' a j j Q'aiQ'Cti: 

S~i = a ] ai \J2S - o) ai a aU 

Sai = S — a* ai a a i, (Bl) 


where a labels sublattice and a = (1,2,3,4, 5, 6 ) and a a i is 
a bosonic operator. Then, keeping only the linear terms and 
using Fourier transformation, the magnetic Hamiltonian can 
be expressed as, 


FIG. 6 : Band structure and DOS of Rb 2 Cr 3 As 3 with experimental 
parameters in the paramagnetic state. 
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H m = ^i?o(l,2)+ffo(2,3)+ffo(3,l)+i?o(4,5) + i7o(5,6) + J ffo(6,4) 
k 

+ 53 Hi (1,4) + H x (1,6) + jfi (2,4) + Fi (2,5) + Hi (3,5) + H x (3,6) 

k 

+ 53 fl2(l, 1) + H 2 { 2,2) + H 2 (3,3) + H 2 (4,4) + H 2 (5,5) + H 2 (6,6), 

k 


Sj SJ\ 3 1 

2) — ~ (cL^CL2k H"~ ®'lkO j 2k) 2 ^2fe^2fc) — tSt/i (Gy£,G.2_fc /c) J^S ’ 

S Jf k 3 

#i(M) = -^ 1 cos(y)(a| fc o 4A: +4 fc ai fc ) - -S J[cos(-^-){a lk a A _ k + a\ k a\_ k ) 

+ S J[(a\ k aik + a| fe a 4 fc) — J[NS 2 , 

-ff 2 (l, 1) = 2SJ 2 cos(& z )(a{ A .a 1 fc) - 25 , J 2 cos(fc z )(a| fc aife + <4 ft au) + J 2 NS 2 . 

In the basis $t (fc) = (a| fc , a\ k ,a\ k , a^ fc , <4 fe , a' 6k , ai_ fe , a 2 _ fc , a 3 _ fc , a 4 _ fc ,a 5 _ fc , a 6 _ fc ), the magnetic Hamiltonian is, 
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A = J X S + 2J[S -2J 2 S + 2J 2 Scosk z , (B8) 

B = \.hS, (B9) 

c = -jj./iS; (bio) 

D = \j[Scos (^), (Bll) 

£ = -|j(Saw(^). (B12) 


The spectrum of the spin wave is given in Fig. 9, with J[ = 
0.7Ji, J 2 = —0.15Ji. In A 2 Cr 3 As 3 , the bond lengths in 
CrlAsl and Cr2As2 planes are different, which leads to dif¬ 
ferent Ji couplings in CrlAsl and Cr2As2 planes. Therefore, 


we plot the spectrum of spin wave (Fig. 4(b)) in the main text 
with Jicri = Ji,JiCr 2 = 0.9J 4 and J[ = 0.7J 4 , J 2 = 
—0.15Ji. Compared with Fig.9, there are gaps on the Bril- 
louin Zone boundary, indicating the intrinsic two sublattices. 
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FIG. 8: Model of the Cr octahedral cluster in A 2 Cr 3 As 3 . 



FIG. 9: The calculated spin-wave dispersion of K 2 Cr 3 As 3 along 
k, direction, with Jicr-i = Jicv 2 = -h , J'\ = 0.7Ji and J 2 = 
-0.15Ji. 






